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Abstract 

In this paper, triple neutral electroweak gauge boson production processes, 
viz. 777 , 77^, 'yZZ and ZZZ productions merged to 1-jet have been studied 
at the leading order in QCD in the context of Randall-Sundrum model at the 
LHC with center of mass energy ^/S = 13 TeV. Decay of Z bosons into lepton- 
pairs has been considered. We present a selection of kinematical distributions 
matched to parton shower and show their deviation from the SM results as a 
result of the RS model. The uncertainties as a result of the factorization and 
renormalization scales are also presented. 
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1 Introduction 


The Standard Model (SM) of particle physics is now well established as the true de¬ 
scription of particles and their interactions at the experimentally accessible energies. 
The recently discovered 125 GeV scalar at the LHC Run-I na, behaves like the 
SM Higgs boson and this hxes the last free parameter of the SM Lagrangian. So far 
there is no indication of exotic searches of beyond standard model (BSM) physics 
and the BSM scales have been pushed further. Remarkable agreement between the 
predicted SM values and the measured cross sections spanning a broad range is a 
validation of the analytical methods and the Monte Carlo tools developed to match 
the challenges on the experimental and theoretical sides. None the less the SM is 
not a complete description of nature on various counts and the issue could only be 
addressed from beyond j3j. 

Run-H at the LHC is now on with higher energies and luminosity. Precision 
measurement of the properties of the newly discovered Higgs boson is a priority and 
will be matched with precise higher order theoretical predictions to look for any 
deviations from the SM. Now that the Higgs boson is discovered it is important to 
look at the massive vector boson scattering (VBS) cross section which is uniterised by 
the Higgs boson. The VBS gets contribution from the non-abelian couplings of the 
electroweak gauge boson sector (a) triple gauge boson coupling (TGC), (b) quartic 
gauge boson coupling (QGC) and in addition, Higgs coupling to the massive gauge 
boson. The TGC will contribute to the di-boson hnal state and the Run-I has already 
placed comparable limits to the anomalous couplings as the LEP experiments. The 
QGC coupling leads to tri-gauge boson hnal states, though the full process would 
involve the TGC, fermion mediated processes and also Higgs mediated processes. 
The QGC has been reported for the hrst time by ATLAS collaboration using the 
VBS process, yielding a hnal state with two same sign W boson in association with 
two jets | 1 ] in a purely electroweak process and also measured the IV 77 production 
cross section [5] which is now accessible with the 8 TeV LHC data set. The CMS 
collaboration has also measured the QGC in the WW'-y and IVZ 7 hnal states [B]. 
So far the observations are consistent with SM predictions and as the sensitivity of 
these measurements improves, the TGC and QGC that lead to tri-boson hnal states 
can not only test the electroweak sector of the SM but also probe new physics. In 
this paper we look at some of the tri-boson production processes merged to 1 -jet in 
the warped extra dimension model. 

In the model proposed by Randall and Sundrum (RS) [7], the non-factorizable 
geometry of the space-time with the inclusion of a single warped extra spatial dimen¬ 
sion, proposes a solution to the hierarchy problem. The SM helds are conhned to 
a 3-brane, whereas the gravity which propagates the full 5-dimensional space-time, 
manifests as massive Kaluza-Klein (KK) modes in 4-dimensional space-time. RS 
model phenomenology of the virtual graviton exchange have extensively been stud¬ 
ied for gauge boson pair production processes viz., 77 [ 8 ], ZZ [9], W^W~ [lO] and 
also for DY production m at the next-to-leading order (NLO) accuracy, because 
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of their rich sensitivity to the model parameters. This in tnrn helps to rednce the 
theoretical nncertainties whereby constraining the RS model parameters. Recently 
the parton shower effect has been considered for those processes in na. Neverthe¬ 
less, the stndy of triple gange boson production processes within this model would 
also be phenomenologically important, as they could effectively participate in inter¬ 
esting new physics searches at the TeV scale. They have been studied in the SM 
at NLO [ini [TH Uni [16] level. NLO results of the triple photon production in the 
SM have recently been presented including the effect of photon fragmentation m 
or matching them with different parton shower (PS) Monte Carlo programmes [TB] . 
These SM processes also serve as potential backgrounds to a number of new physics 
signals coming from different BSM scenarios. For example, the SM 777 process is 
a background to single photon production, together with one techni-pion in tech¬ 
nicolor model, whereas 77 ^’ process in the SM is a background to the signal with 
di-photon plus missing energy in gauge-mediated supersymmetric theories. 

In this analysis, we consider the production of neutral triple electroweak gauge 
bosons in warped extra dimension model at the LHC, i.e., PP VVV X, where 
1/ = 7 , Z and X denotes some hadronic hnal states. Similar processes have been 
analysed at LO in case of large extra dimensional models [19] m m- In fact, study 
of these processes in RS scenario bears equal importance, as their contributions in 
searching new physics using the triple gauge boson productions are undeniable in 
distinguishing physics arising from the potential BSM candidates like supersumme- 
try or technicolor. 

This paper is organised in the following way: in section 2 we present brief de¬ 
scription of the RS model. In section 3 we discuss the merging procedure and 
computational details. Numerical results of various kinematical distributions are 
provided in section 4 and hnally we draw the conclusion of our study in the last 
section. 


2 Neutral Triple Vector Boson Production in RS 
Model 

In the RS model, the non-factorizable geometry is governed by the following 5- 
dimensional warped metric [ 22 ] . 

— rldcj)^ , ( 1 ) 

where 7 ^,^ represents the flat Minkowski metric, (j) denotes the hfth dimension (0 < 
(j) < tt) which is compactihed on a S^/Z 2 orbifold with a radius and k is related 
to the curvature of the AdS^ space-time. Two 3-branes with opposite tensions are 
situated on two hxed points (0 = 0, tt) of the extra dimension. The brane at 0 = 0 
is called the ‘Planck brane’ and the other one at 0 = tt is known as the ‘TeV brane’. 
In the RS scenario, all SM helds are considered to be conhned on the TeV brane, 
whereas the gravity can propagate in the full 4-1-1 dimensions. The interaction 
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among the SM fields with the massive KK excitations of the graviton [231IM] 

is determined by the following Lagrangian, 


-| _ OO 

PI rno 

where cq = is the energy-momentum tensor for the SM 

particles and Mpi is the reduced Planck scale. Note that, the couplings of the 
zeroth KK mode to the SM fields are Mpi suppressed and hence this term can be 
practically neglected. However, the contribution from the higher modes with the 
coupling Co/mo can be of the order of few TeV for a choice of kTc ~ C>(10) [231 ES] 
and they can produce significant observable effects. The masses of the KK mode 
excitations are given by = Xn k, , where Xn indicates to the zeros of the 
Bessel function Ji(x). 

The effective graviton propagator after summing over all the massive KK modes 
except the zeroth one takes the following form m 0. 


where Sij = {pi+PjY, x = y^/ruo and T^ denotes the width of the resonance with 
mass Mn- The total decay width of the graviton can be calculated with the KK 
states decaying to the SM particles [23] [28] in the following way, 

Tn = moclxl_/\n , (4) 

where A„ is given by, 

A„ = Ar + Ay + A*"' + A«" + 5;Ar + 5;A" + Af+ . (5) 

V I q 

Here each corresponds to the coefficient coming from the decay width calculation 
of the process —)■ aa. Unlike the large extra dimension model, the individual 

resonances of the graviton are well-separated in the RS model and they can be 
probed in invariant mass distribution. 

The massive RS graviton could be produced in association with a photon or a Z 
boson and since the RS gravitons also couple to two photons or two Z bosons, the 
SM three neutral gauge boson final state distributions could be altered due to the 
RS contributions and its interference with the SM. 




Sij - MI+ iTnMn 


mn 


— xl) — iXn — 
\ n) "■ mo 


^ - xlf + xl 


( 3 ) 
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3 Merging Matrix Element with Parton Shower 

The leading order neutral triple gauge boson production processes PP —)■ VVV X 
at the LHC come from the subprocess, 

g(pi) + g(P2) V(pa) + V(p4) + V (ps) 

where H = 7 , Z and X is any dual state hadron. This process has been merged with 
the 1-jet process PP —)■ VVVjX in MadGraphS (MG) [29] framework to have a 
better description of different distributions. Due to extra radiation emission, q{q)g 
initiated subprocesses also come up. The merged events are then matched to a 
Parton Shower (PS). The Z bosons are let to decay to lepton pairs, thus accounting 
for the off-shell contributions. 

In LHC, additional jets are often produced from initial state radiation and can 
alter the LO predictions for relevant observables. Generally these additional jets are 
simulated using PS monte carlo. But these QCD radiations in the PS programs are 
generated in the soft and collinear approximation based on Sudakov form factors. 
The widely separated and hard emissions are not well-described in the PS approach, 
whereas the hxed order tree level amplitudes can provide reliable predictions in the 
hard region, but it fails in the collinear and soft limits. Therefore it is also essen¬ 
tial to take into account the tree level amplitude containing additional jets. Both 
descriptions have to be combined in an appropriate matching method by avoiding 
double counting or gaps between samples with different multiplicity. Several algo¬ 
rithms have been proposed for this purpose, mainly based on the event re-weighting 
(eg. CKKW) [30][3T] or event rejection (eg. MLM) [32] . 

The shower-k-r scheme [33], based on event rejection like MLM as implemented 
in MadGraphS is used in this analysis. In this scheme the events are generated 
by MG with a minimum separation in the phase space Qcut and Pr^in between the 
hnal-state partons {ij) and between the hnal-state and initial-state partons {iB) 
respectively which is characterized by the kx jet measure: 

4 = min{pl.,p\)XRl > = p^^ > pl^.^ (6) 

Here = 2 [cosh( 7 j — pj) — cos(0j — 0j)], where 4>i are the transverse 

momentum, pseudo-rapidity and azimuthal angle of the parton i. The kx value 
is set as the renormalization scale at each QGD emission vertex. The events are 
then passed to Pythia[31] for showering. In shower-Zc^ scheme, pythia pr-ordered 
shower is used for showering. Pythia reports the scale of the hardest emission 
(Qhardest) shower and vetoes events based on the kx values of the hardest 

shower emission instead of performing a jet clustering and comparing to the ME. If 
Qhardest > Qcut for lower multiplicity samples, then the event is rejected, whereas 
for highest multiplicity sample an event is rejected if Qhardest > Q^oftesv scale of 
the softest parton in the event from ME. We choose to work with Qcut = PTmin- 

The hxed order merging approach gives a better description of the region of 
hard and well separated jet whereas the parton shower takes care of the infrared 
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region correctly. These merged-matched events provide a realistic framework to be 
compared to the experimental outcomes. The Lagrangian of the RS model is written 
using FeynRules [35] and it is combined together with the SM Lagrangian. The 
universal FeynRules output (UFO) of the combined Lagrangian (he., £ 7^5 + Csm) 
is then imported within MadGraph5 framework and used for the generation of 
events. The model parameters cq and Mi, the mass of the hrst excited KK mode 
have been set as external inputs and we choose to work with the following values: 
Ml = 1.7 TeV and cq = 0.03 which remain within the latest experimental bounds 
provided by ATLAS [SB] |3Z| and CMS [M|[3n] collaborations. In addition, we 
have systematically implemented the KK mode summation algorithm in the spin-2 
HELAS routine [lO]. Kinematical distributions of various observables have been 
recalculated for different di-final states such as, di-photon, Drell-Yan, ZZ, W^W~ 
in fixed order NLO and NLO-I-PS using this present layout and they are found to be 
in excellent agreement with those results that are present in the literature [H [9], [TOl 
mm- This essentially ensures the proper execution of the whole computational 
set-up. We have generated events for the following four neutral triple vector boson 
production processes: (i) 777 , (ii) 77 Z, (m) 7 ZZ and {iv) ZZZ under the above 
mentioned arrangements. The Z bosons are decayed to lepton-pair which will be 
discussed in the next section in detail. Each of these processes consists of three 
types of contributions coming from pure SM, pure RS and the interference between 
these two. 


4 Numerical Result 

In this section, we present numerical results of various kinematical distributions for 
the above four processes. All the results are presented for LHC with center of mass 
energy y/S = 13 TeV. In our analysis, the following set of external parameters are 
used as input: 

mz = 91.188 CeV , sin2(0vv) = 0.222 , 

Gf = 1.16639 ■ 10-® CeV^ , = 132.507 . (7) 

During the generation of events we let the Z bosons to decay to pair. Events 
are generated with loose cuts on the transverse momentum {Pt) and rapidity {y) of 
the final state particles: 

> 15GeU , \y^’^\ < 2.6 . ( 8 ) 

The factorization scale {fip) as well as the renormalization scale (/iij) are chosen as 
the invariant mass of the hnal state particles and MSTW2008LO ( 68 %GL) PDF 
has been used. Throughout this paper, we have considered five massless quark 
flavors (nj = 5) and neglected all top quark contributions. In case of processes 
containing two or three photons in hnal state, a photon separation cut > 0.3 

is used during event generation. AR^^ = {AyY -|- (A 0)2 is the separation of the 
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two photons in the rapidity-azimnthal angle {y, 0) plane. For processes involving 
leptons and photons in hnal state, we have applied > 0.3 and AR^+i- > 0.3. 
For consistently merging 0-jet sample with 1-jet sample we have followed the path 
prescribed in | 33 ] and checked the stability of the cross-section by varying the scale 
Qcut- Additional checks on the smoothness of the distribntions, for example the 
differential jet-rate (DJR) plots, in the ME and PS transition region have been done. 
For the four processes under consideration the choices of Qcut are much the same 
around 90 - 95 GeV where we find best smooth DJR plots and also the matched cross- 
section remains within 13% of the unmatched cross-section. For 777 , 'jZZ and ZZZ 
we choose Qcut = 95 GeV whereas for 77 Z we took Qcut = 90 GeV. Showering is 
done with Pythia P-p-ordered shower as described in the previous section. Different 
analysis cuts used on the final state particles at the time of analysis are described 
in the following subsections for each processes. 

For the processes involving photons, the photons can come from the hard process 
or as a result of fragmentation which is a QED collinear effect. In order to get rid 
of such collinear divergences without involving additional non-perturbative effects, 
the smooth cone isolation criteria on the photons as proposed by Frixione HU , is 

used. A cone of radius R = \J{y — + (0 ~ is considered in the {y — 0) 

plane around a photon satisfying the condition that the total hadronic transverse 
energy E{R) within R < R^ would be less than a maximum limit E{R)max given 
by, 

E{RUax = 6 , P0 f , (9) 

\ 1 — cos R^ J 

where P 0 is the transverse energy of the photon; ey,R^,n are three parameters of 
the Frixione isolation. During event generation we choose = 1,R^ = 0.3, n = 1. 

For the reconstruction of Z bosons from opposite sign lepton-pair, events are 
selected based on the selection criterion: 


\Mi+i- — Mz\ A 15 GcV , (10) 

where Mi+i- is the reconstructed invariant mass of the opposite sign lepton-pair. 
The following transverse momentum and pseudo-rapidity cut on jets during analysis 
in all four processes are used 

F0 > 50 GeV, T]^ < 4.5 (11) 

For all four processes, an extra cut has been put on the final state particles in¬ 
variant mass (M > 600 GeV) for transverse momentum distributions and rapidity 
distributions which are displayed in the respective figures. We have checked that the 
generated events give unbiased results with reference to the choices of generation and 
analysis cuts. Scale dependencies are obtained by varying the renormalization scale 
(pij) and the factorization scale (yp) |12] in the range {jJtR, yp) = (^ijho, ^f/^o), 
where po is the invariant mass of the three vector boson final states or its decay 
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products when Zs are involved. The scale factors kr, kf are in the range (1/2,2), 
we choose the combination {kr, kr) = ( 1 / 2 , 1 / 2 ), ( 1 , 1 ), ( 1 / 2 , 1 ), ( 1 , 1 / 2 ), ( 1 , 2 ), ( 2 , 1 ), 
(2,2) to study the scale variation. The scale uncertainty is represented by taking 
the envelope of all individual variations. 



600 800 1000 1200 1400 1600 1800 2000 2200 2400 


M’’” [ GeV ] 



Figure 1: Invariant masses of 777 (left) and of two hardest photons (right) for 777 
production. 

Next we discuss and present our results for each of the neutral triple boson final 
states separately presenting some select distribution that are of interest for the RS 
model. For all figures, we followed the following convention; we give the distributions 
corresponding to the SM and SM+RS that contribute to the observable for ME+PS 
merged to 1-jet for central scale choice (k^, kr) = (1,1). In the lower inset we put 
the fractional scale uncertainty for ME-I-PS with 1-jet for RS case. 


4.1 777 

Observing the 777 channel has a great advantage over the other triple neutral 
channels, because experimentally it provides a cleaner signature. 777 production 
in the RS model has also been studied in [13]. Here we present the result merged 
with 1 -jet as well as include showering thus improving the result, viable to the 
experimental search. During the analysis level we use more stringent cuts than that 
used at the generation level: (a) cuts on transverse momentum and pseudo-rapidity 
of final state photons used are R/ > 25 GeV, t]"* < 2.5, (b) Frixione parameters 
used are = 0.4, e..y = 1, n = 2 and (c) photon-photon separation cut > 0.4. 
Finally the photons are ordered according to their transverse momentum. 
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Figure 2: Transverse momentum distributions of hardest photon 71 (left) and next 
hard photon 72 (right) for 777 production. 




Figure 3: Rapidity distributions of hardest photon (left) and hardest pair (right) for 
777 production. 


In the Fig. [H we present the invariant mass distribution for tri-photon on the 
left panel and the invariant mass of the hardest di-photon state on the right panel. 
In the invariant mass distribution of the tri-photon the peak appears near 1.7 TeV. 
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The peak is slightly shifted from 1.7 TeV towards higher invariant mass region. 
This is evident from the fact that the RS graviton is produced in association with a 
vector boson. In the expression for invariant mass there are dot products between 
all three momenta of the hnal state photons, thus shifting the peak slightly towards 
the higher invariant mass region. A cleaner signature of RS graviton can be found 
in the invariant mass distribution of the hardest two photon pairs where the RS 
peak appears at 1.7 TeV. From the transverse momentum distributions (Fig. [2]) as 
well as in the rapidity distributions (Fig. |3]), signihcant deviation from SM results is 
observed. In general the 1-jet merged sample gives a harder distributions. For the 
Fig. [T](left panel) the uncertainty at the RS peak is about 7.7% and for the di-photon 
invariant mass (right panel) the uncertainty is about 10%. For the pt distribution 
the uncertainties are larger, for the hardest photon (Fig. [2]) it is about 15% and 
for the second hardest it is about 8.3% around the peaks in the Pt distribution 
which correspond to about half the RS resonance. For the rapidity plots (Fig. [3]) 
the uncertainty is about 10 % in the central rapidity region. 

From the invariant mass distributions (Fig. [T]) an estimation of the signal and 
background events for the LHC Run-II can be made. For this purpose we have 
considered the invariant mass distributions over 1500 GeV where the enhancement 
due to RS graviton signal is clearly visible over the SM background. We hnd that 
for 100 fb~^ luminosity there are 8 signal events over 3 background events from the 
jy /777 distribution (Fig. [T]left panel) and 6 signal events over 2 background events 
from distribution (Fig. [T]right panel). 

4.2 77 ^ 

In 77 Z production the effect of the massive RS KK-modes can be observed in 
the 'j'yZ invariant mass as well as in the di-photon invariant mass (Fig. 0]). The Z 
boson is allowed to decay to pP pair. Minimal analysis level cuts on the hnal state 
leptons and photons transverse momentum > 25 GeV and pseudo-rapidity cut 
< 2.5 are put. In order to fulhll detector resolution, cuts on the separation of 
photon, lepton and jets are imposed: 

R^i>QA, R«>0.4, Rz, >0.7 . (12) 

The invariant mass distribution of 'y'yp~^p~ system shows an enhancement at about 
the RS peak but is additionally shifted compared to the tri-photon case due to the 
mass of the Z boson. In this process the RS graviton can only decay into a photon- 
pair, hence the invariant mass distribution of di-photon pair shows the peak at the 
RS mass. The uncertainty in the enhanced region is about 8 % for the 77 /i■*■/!“ 
invariant mass distribution and for the di-photon invariant mass at the RS peak 
it is about 9.9%. The transverse momentum distributions (Fig. [5)) of each photon 
clearly shows the RS signature appearing nearly at the half of the KK mass. The 
uncertainty for the transverse momentum distribution are much higher, at about 
26%. 
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Figure 4: Invariant mass distributions of 77 /i^/i (left) and 77 (right) for 77 Z pro¬ 
duction. 




Figure 5: Transverse momentum distributions of hardest photon 71 (left) and second 
hard photon 71 (right) for 77 Z production. 


In the Pt distributions in Figs. O [5] the parton shower resums the large logarithms 
in the collinear region which suppress the low Pt cross section. In this region the 
scale uncertainty is also low as the higher logarithms are resummed to all orders. 
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4.3 ^ZZ 
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Figure 6 : Invariant mass distributions of 7 e'*'e (left) and e'*'e (right) 

for 'yZZ production. 

The same set of analysis cuts are used here, as the 77 ^ case. One Z-boson is 
decayed to a e+e“ pair and the other Z-boson to a pair during event genera¬ 

tion. In the 7 ZZ invariant mass distribution (Fig. [H]) there is slight enhancement of 
RS contribution over the SM, around the invariant mass region (1.7 TeV). But the 
effect of massive KK state is best found on the di-Z boson invariant mass distribu¬ 
tion (Fig. [ 6 ]) and it peaks at 1.7 TeV as expected. The pr distributions are not very 
useful to discriminating the RS signatures. The uncertainty for the invariant mass 
distribution of 'ye'^e~fi~ system is about 15% and for the invariant mass of the 
e^e~pZpT system which is as a result of the RS graviton decay is about 8 . 8 %. 


4.4 ZZZ 

For ZZZ production, we choose to decay two Z bosons to two pairs whereas 

the other Z boson is decayed to a e'''e“ pair. For this process also same set of analysis 
cuts has been used, Rf > 25 GeV, 7 ^ < 2.5, where I = e+, e~, p'^, p~. The Z bosons 
are reconstructed according to the criterion described at the beginning of this section 
(see Eq. flTU]) h The reconstructed Z bosons are then ordered according to their 
transverse momentum. Thus in the hgures Zi ^2 represents transverse momentum 
ordered Z bosons. Small enhancement over the SM can be seen in the ZZZ invariant 
mass distribution (Fig. Oleft panel). The uncertainty in that region is found to be 
about 13-14%. In this case also the RS contribution is best found in the hardest 
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Figure 7: Invariant mass distributions of (left) and hardest two 

lepton pairs (right) for ZZZ production. 


two Z boson invariant mass distribution (Fig. [7] right panel). The uncertainty at 
the peak region is about 13.25%. 

5 Conclusion 

In the context of RS model, effects of the exchange of virtual KK graviton have 
been studied on the neutral triple gauge boson production processes at the 13 TeV 
LHC. This process could play a vital role in discriminating physics beyond the SM 
and in estimating the contribution coming from potential BSM scenarios in new 
physics searches. We have incorporated the RS model, using FeynRules in associ¬ 
ation with an algorithm that takes care of the KK mode summation of gravitons 
within MadGraph5 environment and performed a number of checks to ensure their 
proper implementation. 

We have merged P P ^ VVV and P P ^ VVV+j event samples for better pre¬ 
diction of the distributions and observe that it gives harder distributions compared 
to the unmerged sample. To make theoretical prediction closer to the experimental 
situation, we have also matched the merged events with parton shower. Final state 
Z bosons are allowed to decay to either of the following leptonic decay modes: (i) 
Z —>■ (ii) Z —)■ at the time of event generation, thus taking into account 

off-shell effects as well. For process with more than one stable photon, the photons 
are ordered according to their Pt and then the required number of photons are col¬ 
lected, based on their hardness. Likewise, for triple Z-boson production, Z bosons 
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are reconstructed from their daughter particles and then ordered according to the 
hardness. Numerical results of some selective differential distributions for a set of 
kinematical variables have been presented for the merged samples. All these codes 
are flexible enough to incorporate the experimental cuts, different values for model 
parameters etc. and they can be used to obtain numerical results of any kind of 
distributions that would indeed help the experimental collaborations. 

Of the neutral tri-gauge boson hnal states considered here, the tri-photon hnal 
state has the highest rate and can be used to look for signatures of the RS model. 
For the tri-photon process, the invariant mass, Pt distributions of various photons 
ordered in terms of the hardness and rapidity distributions are all good discrimi¬ 
nators of the RS model. The scale uncertainties are by and large within 10% for 
the invariant mass distribution. In the tri-hnal invariant mass distribution the cross 
section is enhanced in the RS resonance peak region, which diminishes in going 
from the 777 to the ZZZ hnal state. The di-invariant mass distribution for all 
four processes is a clear indicator of the RS resonance peak. In the di-invariant 
mass distributions the peak of the RS resonance is most enhanced for the 777 and 
diminishes in going to the ZZZ production process. 
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